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We calculate the Renyi entropy of a positive integer order M for a reduced density matrix of a 
single-level quantum dot connected to left and right leads. We exploit a 2 x 2 modified Keldysh 
Green function matrix obtained by the discrete Fourier transform of a 2 M x 2 M multi-contour 
Keldysh Green function matrix. A moment generating function of self-information is deduced from 
the analytic continuation of M to the complex plane. We calculate the probability distribution of 
self-information and find that, within the Hartree approximation, the on-site Coulomb interaction 
affects rare events and modifies a bound of the probability distribution. A simple equality, from 
which an upper bound of the average, i.e., the entanglement entropy, would be inferred, is presented. 

For noninteracting electrons, the entanglement entropy is expressed with current cumulants of the 
full-counting statistics of electron transport. 

PACS numbers: 05.30.-d, 73.23.-b, 03.67.-a, 72.70.+m 


I. INTRODUCTION 

Full-counting statistics is a powerful theoretical tool 
to investigate the statistical properties of electron trans¬ 
port^. This statistical method enables us to calculate 
the probability distribution of the number of electrons 
transferred between two subsystems, left and right leads 
connected by a quantum conductor. The entanglement 
entropy is also a measure of correlations between the two 
subsystems^—. Suppose we partition our total system 
into complementary subsystems A (the left lead and the 
quantum conductor) and B (the right lead). Then the 
partial trace of a density matrix of the total system p 
over the subsystem B degrees of freedom, 

p A = Tr B p, (1) 

defines the reduced density matrix. The operator of the 
information content, i.e., the self-information associated 
with an outcome described by the reduced density ma¬ 
trix, may be given by I = —In pa (we choose base e). 
The operator is often called the entanglement Hamilto¬ 
nian and its spectrum, the entanglement spectrum^, has 
been widely used to study topological phases. The en¬ 
tanglement entropy is the von Neumann entropy^ of the 
reduced density matrix given as 

(I) = Tr A {paI} , (2) 

where Tr^ means the partial trace over the subsystem 
A degrees of freedom. Technically, it is convenient to 
exploit the “Renyi entropy” of order Af— — . 

S M = Tt a [pa M ] , (3) 

and calculate the entanglement entropy from its deriva¬ 
tive (I) = — liniM->i OSm/OM ( Precisely, Eq. (O is a 
modified definition introduced in Ref. Gl' which is con¬ 
venient for our purpose ). 

The entanglement entropy^ and the Renyi entropy^ are 
closely related to the Levitov-Lesovik formula—, the cur¬ 
rent cumulant generating function of the full-counting 


statistics: They are expressed by a unique quantity, the 
correlation matrix— 4 . However, the relation is limited to 
non-interacting electrons. Recently, Nazarov proposed 
another approach to calculate the Renyi entropy of an 
integer order M by introducing the Keldysh Green func¬ 
tion defined on a multi-contour, which is a sequence of 
M replicas of a standard Keldysh contour—. Ansari and 
Nazarov have further developed this method^— and re¬ 
late the flow of Renyi entropy with the flow of heal—. Al¬ 
though the results are limited to weak coupling between 
the two subsystems, the approach would be promising 
since it enables one to utilize field theory techniques. 

In the present paper, we consider a single-level quan¬ 
tum dot connected to left and right leads [Fig. [Tj and 
calculate the Renyi entropy of a positive integer order 
M by accounting for the tunnel coupling to all orders as 
well as the on-site Coulomb interaction up to the lowest 
order. The reduced density matrix is derived by tracing 
out the degrees of freedom associated with the right lead. 
We will utilize the anti-periodicity of the multi-contour 
Keldysh Green function and perform the discrete Fourier 
transform. The ‘Matsubara frequency ’ 19 : 20 introduced in 
this way 

Xt = n ( 1 ~^w) 1 = (4) 

is a ‘counting field—, which counts an electron transfer 
between replicated Keldysh contours. The resulting 2x2 
modified Keldysh Green function matrix is closely related 
to that previously introduced in the context of the full¬ 
counting statistics [see, e.g. Refs. I2l l2ll - f28l and references 
therein]. This enables us to apply the Keldysh diagram¬ 
matic technique to calculate the Renyi entropy, which 
is formally a ‘Keldysh partition function’ defined on the 
multi-contour. For non-interacting electrons, we relate 
the Renyi entropy with the current cumulants, Eq. (1461) . 
as previously demonstrated by Song et a/— based on the 
correlation matrix. 

Another purpose of the present paper is to examine 
the idea of the full-counting statistics of self-information 
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I. After the analytic continuation of M —y 1 — i£, the 
Renyi entropy turns into the information generating func¬ 
tions^, which is the moment generating function of 
probability distribution of self-information. Although it 
can generate all orders of moments, the statistical prop¬ 
erties of moments other than the first moment m have 
rarely been investigated. We herein calculate the prob¬ 
ability distribution of self-information and find that al¬ 
though the on-site Coulomb interaction weakly affects 
the entanglement entropy, it also affects rare events and 
alters a bound of the probability distribution. We also 
note a simple equality © from which the upper bound 
of the entanglement entropy © can be deduced. 

In the following, we concentrate on the entangle¬ 
ment entropy under the DC source-drain bias volt¬ 
age in the limit of long measurement time. For non¬ 
interacting electrons, our result reproduces the expres¬ 
sion by Beenakker— , Eq. (IfiUl) . In order to estimate the 
accessible entanglement, one has to perform the projec¬ 
tion measurement on electron numbers of subsystems, 
which only generates a sub-leading correction^. 



FIG. 1: The single-level quantum dot coupled to left and 
right leads. We partition the total system into subsystems A 
and B and trace out the degrees of freedom associated with 
the subsystem B. The subsystem A consists of the dot and 
the left lead. The subsystem B is the right lead. 

Here, we re-emphasize two messages, which we feel the 
most important in the present paper; (1) We will intro¬ 
duce the discrete Fourier transform of the multi-contour 
Keldysh Green function defined in Refs. fl5l-|T8l It pro¬ 
vides a feasible way to calculate the Renyi entropy start¬ 
ing from a microscopic model Hamiltonian. It enables 
us to treat electron interactions with a slight extension 
of traditional diagrammatic techniques. (2) Through ex¬ 
plicit calculations, we would like to demonstrate that the 
Renyi entropy may contain the information on the fluc¬ 
tuations beyond the entanglement entropy. We expect 
that the concept of the probability distribution of self¬ 
information could provide a way to interpret the meaning 
of the Renyi entropy. 

The paper is organized as follows. We summarize the 
information generating function in Sec.HTland then intro¬ 
duce our model Hamiltonian in Sec. m In Sec. m after 
we analyze the Renyi entropy for decoupled systems, we 
express the Renyi entropy in the form of the ‘Keldysh 
partition function’ defined on the multi-contour. Then, 
we summarize the discrete Fourier transform of the mod¬ 
ified Keldysh Green function. In Sec. El we present the 
results for noninteracting electrons. In Sec. IVI1 we dis¬ 
cuss the effect of Coulomb interaction within the Hartree 


approximation. Section [VIII summarizes our findings. 


II. INFORMATION GENERATING FUNCTION 


The information generating functio n 29 ! 30 , the moment 
generating function of the self-information /, would be 
obtained from the Renyi entropy © by extending M to 
a complex value 1 — 

Si-it = (,p A (r )“*) = J dIe^P T (I ). (5) 

Here, r is a measurement time. The information 
generating function satisfies the normalization condition 
Si = f dIP T (I) = 1. The n-th cumulant is 


«n> 


d n In S\-j£ 

mr i=o 


( 6 ) 


The first cumulant is the entanglement entropy ©,(/) = 
((/)). The second cumulant (variance) is ((I 2 )) = (I 2 ) — 
(I) . The probability distribution of the self-information 
may be obtained by the inverse Fourier transform, 


PAV = j = (5(1 + In P a (t))). ( 7 ) 


The fluctuations are induced by degrees of freedom as¬ 
sociated with the subsystem B, which have been traced 
out. We note a simple identity, which is reminiscent of 
the Jarzynski equalit y 28 ’ 32 ! 33 , 


(e 1 ) = / d/P(/)e 7 = S 0 = Tr A p A (r) c 


( 8 ) 


which would follow from Eq. © by setting £ = — i. The 
RHS of Eq. © is the rank of Pa(t)° , which is the 
unit matrix in the available many-body Fock space [see 
Ref. 136) . If P(I) > 0, Jensen’s inequality provides the 
upper bound of the entanglement entropy, 


InTrA 


Pa(t)° 


>(I) , 


( 9 ) 


which means that the entanglement entropy cannot ex¬ 
ceed the entropy of the uniform distribution over avail¬ 
able states in the many-body Fock space of subsystem 

A. 


III. MODEL 

The Hamiltonian of the single-level quantum dot con¬ 
nected to left and right leads is 

H = H l + H r + H d + Hr + H int , (10) 

The quantum dot is represented by a localized level with 
the energy e D , 

H d = ^eodlda , ( 11 ) 
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where d a is an annihilation operator of an electron with Tro [pd M ] Tr L [ Pl m ] = n ff SDaMSLaM ■ Spin-resolved 

spin a. The on-site Coulomb interaction is given by, Renyi entropies of the quantum dot and left lead are 


Hint = f/44 d t rf t , (12) 

The left (r = L) and right (r = R) leads are described 
by the free electron gas as 


SDaM =E-/j>( e -P) M ) 

± 

( 20 ) 

SLcrM =IIE^( e Lk) M ■ 

( 21 ) 


k ± 


Hr — ^ ( ^rk^r^rr^rkcr 5 (1H) 

ha 

where a r ka annihilates an electron with wave number k 
and spin a. The tunneling between the dot and the lead 
r is described by 


Hr = E E Jrdlarka + H.c. (14) 

r=L,R ka 


The tunnel coupling broadens the DOS of the quantum 
dot: 


2-it[{u — e£>) 2 + T 2 /4] 


T = T l +T r . (15) 


The coupling strength T r = \J r \ 2 5(uj — e r k) be¬ 

tween the quantum dot and the lead r is assumed to 
be energy independent. The transmission probability 
through the quantum dot is proportional to the DOS as 


T(w) = 2tt 


r L r R 

r 


p(uj). 


(16) 


We assume that initially the dot and the leads are 
decoupled and the Coulomb interaction is switched off. 
Then electrons in each region are equilibrated with the in¬ 
verse temperature /3 = 1/(/cbT) (We set h = fee = e = 1). 
The initial equilibrium density matrix is p eq = plPrPd , 
where 

p D =e -P(HD-y. D £„ 4*0 /Z D , (17) 

p r = e -^ H ^-^^ a ^ a ^)/ Zr , (r = L,R ), (18) 


are equilibrium density matrices of the quantum dot and 
lead r. Here p r (r = L , R. D ) is the chemical potential. 
The equilibrium partition function Z r ensures Tr p r = 1. 
In each region, an electron and a hole obey the following 
distribution functions: 

M = x + e ^- Mr ) ’ /r“ H = 1 - /+ H . (19) 


Let us focus on the quantum dot m and evaluate 
the probability distribution by the analytic continuation 
M —> 1 — it; and then the inverse Fourier transform 0. 
At a finite temperature, we obtain 

PdM) = E 5(1 + hi f±(e D )), (22) 

± 

which satisfies (e I ) = f dlPoM) 6 - 1 = 2 for an arbitrary 
temperature. Although the equality would be valid 
at zero temperature, we have to pay attention to zero 
temperature limits and specify the procedures (I) and 

(II). 

(I) We first take a zero temperature limit for a positive 
integer M = 1,2, • ■ ■ and then perform the analytic con¬ 
tinuation M -P 1 — if;. For the initial equilibrium density 
matrix, the resulting information generating function is 
Si—if = 1 since X)±[/*(w)] M = 0 M + 1 M = 1 (r — L, D) 
at zero temperature for a positive integer M. The prob¬ 
ability distribution is then P T (I ) = 5(1) and the equality 
@ is (e 1 ) = 1. 

(II) We first perform the analytic continuation at 
a finite temperature. For the initial equilibrium den¬ 
sity matrix, since [/^(u;)] 0 = 1 at a finite tempera¬ 
ture, the equality © leads to (e I ) = 2 Na , where 
Na = X)cr I + X)fc<r 1- This provides the maximum en¬ 
tropy of the subsystem A, Na hi 2. The average is the 
thermodynamic entropy, 

(I) =Sa = -J2H fn^D) In f±(e D ) 

G =t 

( 23 ) 

ka ± 

which vanishes at zero temperature. The inequality 0 
ensures that the relative entropy^ between the equilib¬ 
rium distribution and the uniform distribution is non¬ 
negative. 


B. Replica method 


IV. RENYI ENTROPY 

A. Two zero temperature limits 

For a warm-up, let us calculate the Renyi entropy 
of a positive integer order M for the initial equi¬ 
librium density matrix, sm = Tr ld [Tr#p e q A/ ] = 


Following Refs. [l5l - fl8l we formulate the perturbation 
theory of the Renyi entropy 0 of the full-density matrix 
at time r: 

p( T ) = U(r)p eq U (r) f , U(r) = e~ lHT . (24) 

We treat the tunnel Hamiltonian and the on-site 
Coulomb interaction as the perturbation, V = HT + H lnt , 
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and rewrite the Hamiltonian (flUl) as H = Hq + V. The 
Renyi entropy (HJ) includes the partial trace over the sub¬ 
system B , Trs, inside the partial trace over the subsys¬ 
tem A, Tr^. To avoid this complication, we adopt the 
replica method^. We introduce M (M is a positive in¬ 
teger) replicas of subsystem B electron annihilation and 
creation operators: 

a Rka t &Rkam ; a Rka ^ ^ a Rkam ^ ; (^5) 


where the denominator of the RHS is the Renyi entropy 
of the initial equilibrium density matrix. This enables us 
to exploit the Bloch-De Dominicis theorem^ (Appendix 
and the linked cluster theorem, which result in 


i n ^L = y 

SM 


ht 

T7.f 


dt n ■■■dti (TcR(tn)/ • ■ 


n=1 - -cr 

x V(tl)l) M)C 


(30) 


where m = 1, • • ■ , M. Then the Hamiltonian H and the 
density matrix p eq are also replicated by the replacement 
(1251) . We introduce another subscript m to specify m-th 
replicated operators, i.e., H 0m , H m , U m , V rn and p eq m- 
The Renyi entropy (|3|) is expressed by a trace over the 
total system, the subsystem A plus M-replicas of subsys¬ 
tem B as 


Sm =Tr UMPeqMUj^UM-lPeqM-l Um~ i 
X • ■ • X U2Peq 2 U 2 Uip e q 1 U^ 

=Tr 

X • • • X U2lPeq2U2jUijp eq iUi\ 


Um I Peq mUm\Um~ 1 j Peq M —iUm —l\ 


(26) 


where the subscript I indicates the interaction picture 

77 r = p iH OmTTT 
m 1 — ° m • 

The time evolution operator and its Hermite conju¬ 
gate are expanded as U m i = T exp (— i dtV m (t)j) and 
Ujnj = T exp (i fj dtV m (t)i) , where T and T are the 
time-ordering and anti-tinre-ordering operators. Here the 
perturbation Hamiltonian in the interaction picture is 
V m (t)i = e iHomt V m e- iHomt . Then, following Ref.|H we 
introduce the multi-contour C , which is a sequence of M 
replicas of the standard Keldysh contour as depicted in 
Fig. [2] We set a starting point at t = r on the lower 
branch of the first replica C\.— The contour goes to 
Peq 1 at t = 0 along C'i,_ and returns to t = t along 
Ci j+ . Then it connects to t = t on the lower branch of 
the second replica H-- ^ successively repeats until it 
reaches t = t on Cm,+ - Then it connects to the starting 
point t = t on Ci . By introducing the contour ordered 
operator Tc, the Renyi entropy can be expressed as the 
‘Keldysh partition function’, 


Sm =Tr 


Tc exp i j dtV(t)i S jp e q M ---p eq i 

= (tc exp J dtV{t)^j^ 


sm , 


(27) 

(28) 


where the integral over t is performed along the multi¬ 
contour C. In the following, we sometimes write the time 
t defined on C mjS as t ms (s = ±). Then the explicit form 
of the operator V(t)j in Eq. (I?51l is V(t ms )i = V m (t)i. 
Note that the contour-ordering operator Tq also acts on 
Peqm residing at t = 0 m ±. The normalized expectation 
value is defined as 


(29) 


The subscript c means that only connected diagrams are 
taken into account. Equation (1301) is the starting point 
of the following calculations. 



FIG. 2: Multi-contour C, which represents successive con¬ 
nections of the M replicated Keldysh contours, starting at 
t = t on the lower branch of the first replica C 1 ,-. Shaded 
boxes are M replicas of the initial equilibrium density matrix 
Peqm (m = I, - - , M ). The solid line connecting C 2 ,- and 
Cm- 1 ,+ indicates a ‘greater’ Green function + ’ 2 ~ (t, t'). 


C. Modified Keldysh Green function 


The diagrammatic expansion of the Keldysh partition 
function (1281) is performed based on the multi-contour 
Keldysh Green function^— . We relegate the details to 
Appendix [B] and summarize the multi-contour Keldysh 
Green function for an electron in the left lead, which is 
a part of the subsystem A. This is a correlation function 
of a Lka 011 C m',s' and a L kcr on C m y. 


/, ,/ \ ms,m s /, ,/\ 

gLka\tmsit m 's' ) ~9 Lka Fd ) 

= — i (Tc a Lka(tms)I^Lka^m' s')/ M ' 

(31) 


This is a component of a 2 M x 2 M Keldysh Green func¬ 
tion matrix g Lka(t,t') [See the explicit form Eq. (1B3I) in 


(0)m — Tl [&PeqM ' * ’Peql] /sm j 
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Appendix [B]. By exploiting anti-periodicity in the repli¬ 
cated Keldysh space, we perform a discrete Fourier trans¬ 
for m 19 ’ 20 , 

M -1 

s x L w)= Y, fe Lk a {t,t')] rntrn ,j*^ m - m, '>. 

m—m'—O 

(32) 


The inverse discrete Fourier transform is 


M—l 

fe =jj E 

e=o 

(34) 


This is 2 x 2 Green function matrix defined in a single 
Keldysh space as 


g Lk<r(t,t') 




(33) 


The parameter Af ©, the ‘Matsubara frequency ’ 19 ’ 20 , 
satisfies Am-i-< = —\e- This parameter is the count¬ 
ing field for electron transfer between different replicated 
Keldysh contours. The explicit form of Eq. (1551) is 


g X Lk At,t') = 




(t-f) - ft\(eLk)0(t' - t) 
-fE.\{£Lk)e~ iX 


ftx(tLk)e iX 


where the modified electron and hole distribution func¬ 
tions are 


/L,A( £ )= 1 + e /3 (e -^)+*A , /l,a( £ ) = 1 - /E( £ ) • ( 36 ) 


In the limit of zero temperature, they are independent 
of the counting field A, lim^oo A (e) = 6(±(hl — e)). 
Then Eq. (155)) becomes the ‘modified Keldysh Green 
function’ introduced in the theory of full-counting statis- 
tics^i^Ir— . This fact enables us to relate the Renyi en¬ 
tropy with the full-counting statistics in a novel manner, 


which does not rely on the correlation matrix — . The 
bare modified Keldysh Green function of the quantum 
dot is given in the same way [see Eq. (IC3I) in Appendix 

EJ. 

For an electron in the subsystem B , the right lead, 
replicated annihilation and creation operators, aRkam 
and o/ Rk , reside only on the same m-th Keldysh con¬ 
tour C rn ,± [see Fig. [5] The contour C m ,± starts at 
t = T m _, goes to t = 0 m _ = 0 m>+ along C m _ and 
returns to t = T mi + along C m ,+\- The multi-contour 
Keldysh Green function is 


J 


SST * (M') = -i {TcankAtmsW Rka (t' ms ,) I ) M = -*Tt 


Cm O'Rkcrm ms )l a ^Rk crm (t'ms')lPR ,1 


J m,m • 


(37) 


Here, Tc m is the time-ordering operator along the contour C m . The 2x2 sub-matrix of the 2 M x 2 M Keldysh Green 
function matrix is [gRkcr(t, = g Rk^{t,t')S mim >, where 


&Rkcr {t-j t ) — 




f R {tRk)0(t - t') - /+ {e Rk )6(t' - t) 
~ fR i e Rk) 


Ir (tRk) 

Ir {tRk)0{t' ~ t) - /+ ( e Rk )d(t - t') 


(38) 


r 


This is the same as the standard Keldysh Green func- V. NONINTERACTING ELECTRONS 

tion except for the minus signs at off-diagonal compo¬ 
nents, which are attributed to a different choice of start¬ 
ing points 

A. Linked cluster expansion 


The linked cluster expansion Eq. (EH) for the non¬ 
interacting case U = 0 can be done straightfor- 
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FIG. 3: The solid line indicates a ‘greater’ Green function of 
a replicated fermion s ( t , t'). Since it connects the same 

replica (m — m/ — 1 ), by deforming the closed time-path, it 
is reduced to the standard Keldysh Green function except for 
the choice of a different starting point for the Keldysh contour. 

wardl y 21 ! 23,38 : 

In — = V l n Ei t (39) 

SM a S a M 

In = -Tr [g D a^a] - ;:Tr [g Da T, a g D a^cr} -■ 

SaM l 

(40) 

The self-energy appears after tracing out the degrees 
of freedom associated with the left and right leads: 
£ a(t,t') = Yjrk\M 2 9rka(t,t'). The product and the 
trace in Eq. (TTO1) should be understood as the integral 
along the multi-contour C, e.g., 

Tr [g D a^a] = f dt 2 dtig Da (ti,t 2 )E a (t 2 ,ti ). (41) 

JC 

After we project the time defined on C onto the real axis 
and perform the discrete Fourier transform, we obtain 

o m-i 

In— = E Wa T (\t), (42) 

S ° M 7^o 

War(X) =Tr ln [l - g^E^] , (43) 

where 1 = diag(l,l) is a unit matrix and T 3 = 

diag(l, —1) is a Pauli matrix in the 2x2 Keldysh space. 
The trace is performed over the real time t € [0, r] and 
the 2x2 Keldysh space. The modified self-energy is 

E = \Jr \ 2 E gRka + \ Jl I 2 E ^Lka ■ 

k k 


At zero temperature, Eq. (1431) becomes the current cumu- 
lant generating function of non-interacting electrons (see, 
e.g. Refs. [2lll23ll24ll28l and references therein). Therefore, 
Eq. (1T21) connects the Renyi entropy and the full-counting 
statistics. 

In the remainder of this section, we consider the zero 
temperature limit (I) in Sec. II V Al The current cumulant 
generating function is expanded in powers of iX as 

(45) 

o n= 1 

where C Tj „ is a n-th current cumulant. By plugging 
Eq. <1351 into Eqs. (13U1) and fl^l) and by using the re¬ 
lation A” = 0 for odd n, we obtain 



(46) 

This equation is consistent with the results of Song et al. 
[Eqs. (2.24) and (2.25) in the supplemental material of 
Ref. Q. Further calculations lead to 

00 n /o \ 2 n 

1 ^M=E(gj(-ir b (C(-2n, ( 1 -M)/2) 

- C(-2n, (1 + M)/2)), (47) 

where £(s, a) = )Ty o ( a 4" ^) _s is the Hurwitz zeta func¬ 
tion. The entanglement entropy obtained from this equa¬ 
tion formally reproduces the results of Klich and Levi- 
tovi, as demonstrated in Ref. H, 

00 _ 

(I) = E 7 ^r(-l) T,+1 (27r) 2 "B 2n , (48) 

71=1 ' 

where B n is the Bernoulli number. 

Let us consider the Gaussian approximation, i.e., we 
keep only the lowest cumulant, the second cumulant C Tj2 . 

n M- l 2/1 \ 

In S m = ^T,^ 2 = C ^y(m- M ) ■ ( 49 ) 

e=o ' ' 

The entanglement entropy is (I) = C Ty2 ir 2 /3. The proba¬ 
bility distribution of self-information I obtained from the 
inverse Fourier transform 0 is 

P T (J) = e -< / > / 2 <5(/ - (/)/2) - 

\J2I / (I) - 1 

x Jr (i(I)/2y/2I/{I) - l) , (50) 

where J n {x) is the Bessel function. The self-information 
is almost exponentially distributed and the lower bound 
is the half of entanglement entropy / > (I) / 2 . 

Recall that we are interested in the situation when the 
source-drain bias voltage is applied and thus the average 


( 44 ) 
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grows linearly in r, (I) oc r. In this situation the inverse 
Fourier transform 0 can be done within the saddle-point 
approximation, i.e., the Legendre-Fenchel transform^, 

In P T {I) ss min (In —i£*I) , (51) 

where £* is a pure imaginary number. This results in 

lnP T (/)»(/>V2//(/)-l-/. (52) 

Figures [4] (a) and (b) are the information generating 
function on the imaginary axis and the probability dis¬ 
tribution (15211 . respectively. The information generat¬ 
ing function is defined on the domain —oo < < 1 

[panel (a)]. The most probable value, the peak posi¬ 
tion, is given by the average, i.e., the entanglement en¬ 
tropy [panel (b)]. The duality property of Legendre- 
Fenchel transform relates the slope of the logarithm of 
information generating function I* = d,^ In (prob¬ 

ability distribution i £* = —dj In P T ) with the argument 
of the probability distribution In P T (I*) (the information 
generating function lnS'i-jg.)^. The logarithm of in¬ 
formation generating function is convex and behaves as 
ln5i_j£ ss ((/)/2)i£ for *£ —» —oo. This indicates that 
the lower bound is (I)/ 2. It diverges at the boundary 
= 1 [panel (a)]. The divergence, in turn, shows that 
the large fluctuations are not bounded and In P T {I) ~ —7 
for / —»• oo. The divergence at = 1 implies that the 
equality © is not well defined and we have to go beyond 
the Gaussian approximation, as we will see in the next 
section. 


B. Limit of long measurement-time 

In the limit of long measurement time r —>• oo, the ex¬ 
tensive component of W aTl i.e., the scaled current cumu- 
lant generating function, Jg (T (A) = linv^oo W crr (A)/r, 
is relevant. In this limit, Eq. (l43l) is calculated as (Ap¬ 
pendix [D]) 


Fgo{ A) =-^ J duj In , 

(53) 

ocq f£(u) + fE(u>)e^ 

(54) 



FIG. 4: (a) The information generating function on the 

imaginary axis within the Gaussian approximation. The 
information generating function is defined on the domain 
—oo < < 1. (b) The corresponding probability distri¬ 
butions of self-information for (7) 1. 

generating function of the binomial distribution with an 
energy dependent transmission probability: 

1 r^ L 

A) = — / duj In [1 + 7»(e lA - 1)] . (56) 

Zn “Vi? 

Here we consider the positive bias voltage n = hl~ Hr > 

0 . 

Using Eq. (15^1) . the spin-resolved Renyi entropy is cal¬ 
culated as (Appendix [Dll 


where a trivial constant is subtracted in order to satisfy 
the normalization condition J-Go-(O) = 0. The effective 
electron and hole distribution functions, and fj = 
1 — f~li are bounded to the interval [0,1]. The former is 
the effective transparency from the right lead to the left 
lead: 


/t (w) = Tfw)^(w) + 7J(w)/^(w), (55) 

where the reflection probability is 1Z(uj) = 1 — T(w). At 
zero temperature, we obtain the scaled current cumulant 


„ M—l 

In —ZM. F g <t{\ t) 

S ° M 7^o 


=T 


dco ( fj{w) M + fZ{w) M \ 
2 - + 


(57) 


In the limit of zero temperature (I) in Sec. IIV Al it be¬ 
comes 


rtr. i 

In S aM ~r — In [T( co) M + TZ(uj) m ] . (58) 

•Vr ZlT 















In the extended wide-band limit=£, where the level broad¬ 
ening is large enough, fx <C T, or the dot level is far away 
from the Fermi energy, \ejj — fi r \ T, p,, the transmission 
probability is energy independent: T(w) ~ T(0) = T- 
Then the Renyi entropy is as follows: 

s M = n s - m ~ ( rM + nM ) JVatt - ( 59 ) 

G 

where iVatt = 2r/r/h is the number of attempts, i.e., the 
number of electrons injected into the quantum dot from 
the left lead during the measurement time r. When N att 
is a positive integer, Eq. is the relative information 
generating function of the binomial distribution^. The 
entanglement entropy 

(I) = -N att {TlnT +Kin'll) , (60) 

reproduces Ref. From Eq. ©> we obtain the num¬ 
ber of available states in the subsystem A , (e 7 ) = 2 Nm . 
The available states are limited to the Fermi window 
Hr < oj < hl, since at zero temperature, the electron 
states outside this window are empty or occupied. The 
inequality © becomes 

JV a tt(ln2 + TlnT + ftlnft) >0. (61) 

For the positive integer N a tt , the inverse Fourier trans¬ 
form of Eq. (15U)) can be done analytically: 


N a, 


Pr(D = E 


A U! 


n —0 


n\(N att - n)\ 


r n iz 


N att Tl 


x 5(1 + n In T + (iVatt — n) In 1Z ). (62) 


For the binomial process, the distribution is symmetric 
when we switch the transmission probability and the re¬ 
flection probability T TZ. Figure [5] (a) shows the in¬ 
formation generating function on the imaginary axis for 
various T(< 1Z). They satisfy the normalization con¬ 
dition In Si = 0 and the equality © In So = IV att In 2 
independent of the transmission probability. Figured] (b) 
shows the corresponding probability distributions calcu¬ 
lated within the Legendre-Fenchel transform (15B 22. In 
the limits of if; —> oo and —>• —oo, the information 

generating function behaves as In Si_,{ ~ (~N at t lnT) 
and lnSi_,{ « (—IVatt In72.) respectively [Fig. [5] (a)]. 
Therefore, the lower and upper bounds are I min = 
—IV a tt In 1Z and / max = —iV at tlnT, which is consistent 
with the analytic expression (1621) . The upper (lower) 
bound corresponds to the self-information of a sequence 
of ./Vatt events where all _/V a tt injected electrons are trans¬ 
mitted (reflected). The probabilities to find I = I m i n and 
I = fmax are given by P T (/ min ) = 1Z Nm = exp(-/ min ) 
and P r (/ max ) = T Nm = exp (-/ max ), respectively. Fig¬ 
ure [5] (c) shows the lowest 4 cumulants as a function of 
the transmission probability T. The higher cumulants 
((/”)) (n > 2) increase around T ~ 0 or 1 and vanish at 
T = 1/2. The distribution takes a simple form, the delta 
distribution, P T (I) = 6(1 — N att hi2), at T = 1/2. 



<5 



FIG. 5: (a) Information generating functions on the imag¬ 

inary axis and (b) probability distributions for T = 0.1, 0.3 
and 0.495. The upper and lower bounds are on lnP T = —I 
(dotted line), (c) The lowest 4 cumulants as functions of the 
transmission probability. 


Let us go back to Eq. m and calculate the average 
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at a finite temperature. 

i 1 ) = ~ ^ j d ^Y ft M ln ft M 

+ s a + t^ J (63) 

The first line is the contribution from electrons fluctu¬ 
ating at the boundary between the subsystems A and 
B. The second line is the thermodynamic entropy mi¬ 
nus the over-counting term. Since we can derive Sq = sq 
from Eq. m, we check that the equality © ensures 
(e I ) = 2 NA . 

VI. COULOMB INTERACTION 


The Hartree diagram is depicted in Fig. |6](c). We assign 
two counting fields Xg and Xg' to two loops and perform 
summations over both of them. 



<7 A <t A o" A o" A 

(b) —*— = —*— + —«■—(D—*■ 

aX. 



Here we calculate the correction induced by the 
Coulomb interaction by exploiting the Keldysh diagram 
technique. For non-interacting electrons, the series ex¬ 
pansion, Eq. (17(71) . corresponds to the diagrams depicted 
in Fig. [6] (a). Each thin solid line represents the bare 
modified Keldysh Green function and each circle 
represents the self-energy For non-interacting elec¬ 
trons, the diagrams consist of a single closed electron 
loop and thus only a single discretized counting field Xg 
appears. In the presence of the Coulomb interaction, we 
have to account for diagrams consisting of more than two 
electron loops, which carry different discretized counting 
fields. Consequently, the link between the Renyi entropy 
and the full-counting statistics Eq. (1771) does not hold, as 
we will demonstrate in the following. 

Let us calculate the on-site Coulomb interaction cor¬ 
rection up to the lowest order in U. The first order ex¬ 
pansion in U is 

g(!) r 

In ~iU / dtg D -f{t,t)gDi{t,t) 
sm J c 

M r r 

=«EE S 

m=ls=± 

■TJ M 1 pT 

=J7 Y Y s dt 9ni aB (t,t)9Di ,aa (*>*)■ 

M er=o s=± 

(64) 


By replacing the bare modified Keldysh Green func¬ 
tion g/,^ with the full modified Keldysh Green function, 
which follows from the matrix Dyson equation [Fig. [6] 

(b)], 

CL = SDc r + SDa T ^a T 3^Dr7 > (65) 

(Appendix [C]) we obtain the Hartree term, 


c(l) 
ln Am_ 
sm 


iU_ 

M 


M -1 


E S / dt *) G Di SS ^ *) • 

</=0s=± J ° 



FIG. 6: (a) Diagrams for the logarithm of the “Renyi en¬ 

tropy” m of non-interacting electrons Eq. pOl). The thin 
solid lines represent the bare modified Keldysh Green Func¬ 
tion gEq. (IC3I) . Circles represent the self-energy E**, 
Eq. ©, which is related to the degrees of freedom of the 
reservoirs. The discretized counting field A^ is conserved on a 
single electron loop, (b) Dyson equation for the full modified 
Keldysh Green function G^, Eq. (1651) . (c) Hartree diagram. 
Two loops carry different discretized counting fields, Xe and 
Xg>. 


The Hartree term in the limit of long measurement 
time reads 


cf 1 ) 

In—— «tMV' U5n 9M n a m, q, (67) 

71/T ZJ 


where a =7 (),) for a =| (t). The classical component 
Sn S M and the quantum component n a M,q of dot electron 
occupancy are given by 




2 iM 

e=o 

r fl / a 1/2 H - Sn&M ? 


£—0 

7 r\ i ^(7 AT 

= H-^D ln - 

Mt S aM 


M 


( 68 ) 


(69) 


The Hartree term is interpreted as a correction caused 
by M-dependent renormalization of the dot level, eu —>• 
cd + USnsM- For M = 1, the quantum component van¬ 
ishes n a i tQ = 0, as we see from Eq. 153). The classical 
component becomes 5n s i = rig — 1/2. The spin-resolved 
occupancy is as follows [hereafter, we concentrate on the 
zero temperature limit (I) in Sec. IIV A] : 


Tlfj 


1 

2 


+ E 



fl r €]J 

r/2 


(70) 


(66) 
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The M-dependent correction is 




r L -T R 

2T 


r^L 


duip(u>) 


'MR 


TZ(uj) 


M -1 


T(w) M +7 Z(u) 


M 


-1 


(71) 


It turns out that the Hartree term affects rare events, 
which correspond to it —> ±oo limits after the analytic 
continuation M — > 1 — it- When the condition T(oj) > 
1Z(uj) is satisfied in the Fermi window pr < oo < pp, the 
classical component in each limit reads 


lim 5n s ,i-i(, 

i£—}—oo 


lim Srias-it 

i £—»oo 


= E 


i 

— tan 
2ir 


p r — 6d 

r/2 


5n7 m ■ 


(72) 


= E 


— tan 
27T 


p r — CD 

(r r - r F )/2 


r sym 

On/ . 


(73) 


where f = L(R) for r = R(L). The former is the dot 
occupancy (subtracted by 1/2) for the symmetric cou¬ 
pling r L = r^. The classical component Sn&p-ig can 
be simplified further in the case that the dot level is be¬ 
tween two chemical potentials, pr < cd < Pl, near the 
symmetric coupling, |r L - I\r|/2 <C \p r - £d\ -C T/2: 

f 0 (it —)> — oo) 

Sn^i^ » { 0 (i£ = 0) . (74) 

[ sgn(r L - T r )/2 (it; —>• oo) 

This indicates that for the upper bound i£ —>• oo, where 
all injected electrons are reflected, each electron observes 
that the dot is fully occupied for or empty 

for Tfl > r^. In contrast, for the lower bound i£ —> 
—oo, where all electrons are transmitted, each electron 
observes that the dot is half occupied. 

Figure [7] (a) shows the £ dependence of the dot occu¬ 
pancy 5n a i-i£ for the nearly symmetric coupling k 
T/j. At U ; —» —oo, where all electrons are transmitted, 
the dot occupancy is slightly modified as compared with 
that at i£ = 0. On the other hand, at it; —>• oo, where all 
electrons are reflected, the dot occupancy is enhanced 
(r L > r fl ) or suppressed (r L < r fl ). Figure [7] (b) 
shows the corresponding probability distributions. The 
upper bound is modified as compared with that with¬ 
out Coulomb interaction (dot-dashed line) although the 
peak position, i.e., the entanglement entropy, is almost 
unchanged. 


VII. SUMMARY 

In summary, we studied statistical properties of infor¬ 
mation content in the presence of a Coulomb interaction. 
We calculated the Renyi entropy of a positive integer 
order M by exploiting the multi-contour Keldysh Green 
function. We demonstrated that at zero temperature, the 
discrete Fourier transform of the multi-contour Keldysh 






FIG. 7: (a) Dot electron occupancies as functions of it 

and (b) the corresponding probability distributions for r l = 
0.5ir (solid line) and r r = 0.49r (dashed line). The dot 
level is off the symmetric point pl — £d ^ £d — Pr- The limit 
it —>• oo (it —>• — oo) corresponds to a rare event in which all 
electrons are reflected (transmitted). The dot-dashed line in 
panel (b) is the result without Coulomb interaction [7 = 0. 
Parameters: pl = ~Pr = —tD = 0.5T, U = OTT. 


Green function is compatible with the modified Keldysh 
Green function introduced previously in the context of 
full-counting statistics. For non-interacting electrons, we 
relate the current cumulant generating function of the 
full-counting statistics, the Renyi entropy and the en¬ 
tanglement entropy without relying on the correlation 
matrix. We further calculate the probability distribu¬ 
tion of self-information by the inverse Fourier transform 
of the information generating function obtained by the 
analytic continuation M —>• 1 — it- Within the Hartree 
approximation, we demonstrate that, in the vicinity of 
the perfect transmission, the dot occupancy is modified 
for rare events. Consequently, the upper bound of the 
probability distribution of self-information is modified. 
We point out the equality reminiscent to the Jarzynski 
equality, from which the upper bound of the entangle¬ 
ment entropy could be obtained. 

In short, we feel there are two important and concrete 
messages in the present paper; (1) The discrete Fourier 
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transform of the multi-contour Keldysh Green function 
in the replicated Keldysh space provides a feasible way 
to calculate the Renyi entropy starting from a micro¬ 
scopic model Hamiltonian. It enables us to deal with 
electron interactions by a slight extension of traditional 
diagrammatic techniques. (2) The Renyi entropy may 
contain the information on the fluctuations of informa¬ 
tion content beyond the average value, the entanglement 
entropy. The concept of the probability distribution of 


self-information provides one way to interpret the mean¬ 
ing of the Renyi entropy. It could be interesting to 
consider the conditional joint probability distribution of 
energy — and information content and analyze the pos¬ 
sible connection between these quantities^!. 

We thank Dmitry Golubev, Ryuichi Shindou and 
Kazutaka Takahashi for their valuable input. This work 
was supported by JSPS KAKENHI grants (grants no. 
26400390 and no. 26220711). 
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Appendix A: Bloch-De Dominicis theorem for the normalized expectation value 


We demonstrate the Bloch-De Dominicis theorem!! for the normalized expectation value (EHll through calculations 
of the following 4-th order term as an example (in this section, we neglect the subscripts k and er for simplicity). 

(-i) 4 ^a L (f 4 )}d(f 4 )/d(f 3 )|aL(t3)Mi2)la fl (f 2 )/a fl (fi)}d(fi)/^ , (Al) 

where t 4 > f 2 £ Ci,_), t 3 £ C m0j+ , and f 4 £ C mo+mi+ [Fig. [5]. This is calculated as 


Tr 


( RlPd) M mi m °a^ L d(pLPD) mi d^aL(pLPD) m ° Pmd^a R ia) R1 d 


e -i6 D (t4-t 3 -t2+ti)-iei,(t3-t4)-ie R (t 2 -ti )/ SlmSdm ) 

(A2) 


The density matrices of the replicated subsystem B , p Rrn disappear except for p R \. Following the standard proce¬ 
dure^, we calculate the trace in Eq. (IA2I) as 


Tr p D M mi ^dpD^d^pD m °pRid^aR\a} R1 d fL,M-m 1 {^L)sLM 
=Tr [p D M - m '- m °dp D m 'd)p D m °dU] fn{e R )f L , M -mAe L )8LM 

- [fD,m i ~ fD,M-m 0 (efl)/t),m 0 +mi ( e £>)] Ir ( e R)fL,M-m x (e l)slmSdM ■ 


(A3) 


where the modified Fermi distribution function is defined in Eq. (IB2I) . Then Eq. ED is expressed with multi-contour 
Keldysh Green functions (IB3I) as 

-^-' mo+ (t 1 ,t 3 )ffr + ’ mo+mi+1 -(i3,t4)ffr +mi+1 -' 1 ”(i4,t2)<7fl-’ 1 “(i2,G), (A4) 

which proves the Bloch-De Dominicis theorem. The first and second terms correspond to the diagrams depicted in 
Fig. [8] (a) and (b), respectively. 


Appendix B: Discrete Fourier transform 


In this section, we summarize the multi-contour Keldysh Green function introduced in Ref. 13 Let us calculate 
Eq. (13111 for s = +, s = — and m > m' as an example. Paying attention that the replicated equilibrium density 
matrices p eq m (m = 1, • • • , M) also obey the contour ordering operator Tc, Eq. (T5TT) is calculated as 


dLkc T (t, t ) — *Tr TcaLko-(tms)l a \ J ka^'m l s')lPeqM ' ' ' Peq 1 / SM 


= — * Tr 


p L M- m+m '- 1 aLk a(t) I p L mr “ m + 1 a[.( t , ) I \ /s LM = -ifL^m-m'+licL^e-^-^ , (Bl) 


where slm = Tier s LaM■ Here, the Fermi distribution function is modified and M-dependent; 

e -mp(e-Hr) 

1 4 . g -MP(e-n r ) 


fr,m (^) 


(B2) 













FIG. 8: Diagrams of the 4-tli order term, Eq. ED- 


We calculate the other components in the same way and obtain a 2 x 2 sub-matrix of a 2 M x 2 M multi-contour 
Keldysh Green function matrix connecting branches C' 7Tlj ± and C m i ± as 


[g Lka(t, t )] m>m / 


9 L ka \t,t) 9 L ka (*>*) = A e -i€Lk( 

SiW (M) (M)J 

h ,771—771' ( e Lk) 
fL,m—m' — 1 (e_Lfc) 

fLfi{tLk)0{t - t') - fL,M^ e Lk)0{t' — t) 
—Jx,M-i(eLfc) 

~Sl ( e Lk) 

f —m' — l (e L k) 


-O 

f L,m— m'+l i^-Lk) 
fL,m—m' ( e Lk) 

/i,i(eife) 

fL,o{ e Lk)0{t' ~t)- fL,M( e Lk)d{t - f') 
f-\- 1 (e Lk ) 

~fL,M+ ra— m' (eifc) 


(m > to') 

( to ' = to ) 

(to < to') 
(B3) 


The multi-contour Keldysh Green function (11 hill pos¬ 
sesses the discrete translational symmetry: 

[gLfccr] m _ m / [gL/cer] m;m / 5 (B4) 

where to — m' = —M + 1, • ■ - , M — 1. It also satisfies the 
anti-periodic boundary condition: 

[gLfccr] n = — [gLfeo-] n + JVf ! (B5) 

where n = —M + 1, • • • , — 1. Therefore, it is expedient 
to set [gLfc<r] ±M = — [gifccrJo and continue the domain of 
the function [gLfc CT ] ra periodically with the period 2M-. 
The discrete Fourier transform and the inverse discrete 


Fourier transform read 

1 2M—1 

£ [gLfe<r] n e i27rjri// ( 2M ), (B6) 

71=0 

1 2M-1 

[g^]„=^ £ g ife ab>- i27rj ' n/(2M) . (B7) 

j=0 

By using the anti-periodic boundary condition (IB5I) . we 
rewrite Eq. (1B6I) and demonstrate that an even j compo¬ 
nent vanishes: 

i i I /_i \j 

g Lk,\j] = £ £ j [g LkAn e i7tjn/M , (B8) 

71=0 

Then, by setting g£ = g Lka [M +1] {l = 0, • ■ ■ , M - 1), 
Eqs. (IB6I) and (IB7I) are reduced to Eqs. (l3lil) and (IMll . 
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In order to calculate the discrete Fourier transform where X(. = e l7T ( 2( ' : + 1 )/ M , The discrete Fourier transform 
(1321), we fix m! and rewrite the summation as of g^ m '~ is calculated as 

M 

SlUM 7 ) = E , (B9) 

m =1 


j 






m - ■ ■ ■ — fL,M x i 1 + /i, 1 + fL,2-m' x l + -b m ) 

M—l 

E fL,j+n4 = -ie- iek(t - t ' ) fl iXi (e Lk )e lXe . (BIO) 

3= 0 


The causal component is calculated as 


+ (*I *0 = - ie J£fe(t ‘ } m - fL,M-l x e 1 + /l,0^(* - t') - f ~ t) + f L , 1%\ + 




M—l M—l 

0(t - ?) E fL.M + W - t) E hj+iXe 

3=0 3=0 


.3 +1 


_ _ i e 


/i ,~ *0 “ /l, \ e ( € k)0(t' - t) . 

For the other two components of Eq. (1551) . we repeat the same calculations and obtain Eq. (1551) . 


(Bll) 


Appendix C: Full modified Keldysh Green function for the quantum dot 

Here, we solve the matrix Dyson equation (1651) in the limit of t —> oo. The Fourier transform of the bare modified 
Keldysh Green function of the left lead (1351) is 


gEM = / d{t-t')e iu(t 4 ' } g LkAt,t’) = 




-2Trif+ x (uj)e lX 6(u} - e L k) 

trttZMe**!* -«») -SSl - 


oj+ir;—e_n, ui—ir]—eLk 

?L, \( 


ui—ir)—£Lk u+iv—£Lk 


(Cl) 


where 77 is a positive infinitesimal and the delta function is defined as S(ui) = r]/[n(oj 2 + r/ 2 )]. The modified self¬ 
energy dMl) is21— , 


siM = -i £ T 


r=L,R 


1-2EH 2/+ A >)e : 


iA r 


r,A r 


-2/rTA.He-^ 1-2/EM 


, (Al = A, A_r = 0). 


The bare modified Keldysh Green function of the quantum dot is given in a similar form as Eq. (El: 
g x D<T (t,t') = 


f D ,x( e D)0(t - t') - f£ tX {e D )0(t ~ t) /EM)e a 

-/w(fD)e _ii /p.aMM*' -1 ) - fnAe D )0{t - t') 


The full modified Keldysh Green function obtained after solving the matrix Dyson equation (1551) is 


GdctM — 


D a (w) 


-1 


(w — ep ) 2 + r 2 /4 _ * Er Br/ r _\ r (w)e 


- iXr ' ep - w - i E r r"[l/2 - /+ M] 


-ep-*E r r r [l/2 - /+ (w)] -i Er r r/EMe 


=1 + 7»[/+ A /* (« <A - 1) + f+fl A (e"“ - 1)] = 

/L( w ) + /i( w ) e 


(C2) 


(C3) 

(C4) 

(C5) 


Appendix D: Details of calculations in Sec. IV 51 TrlnG^ CT 1 gp cr - By using Eqs. (IC1I) and (IC4I) . the 

By using the full modified Keldysh Green function, 

Eq. (1651) . Eq. (1451) can be rewritten as Wen-(A) = 
























14 


limit of t — > oo is calculated as 


W (7T (A) / dw In 


27T 


detG^ CT (w) 

.-l 


det G 


A—01 
Da ' 


W 


S 


cku In 


dct G^-°(o;) 

4 -l 


detg^ CT (w)" 

= ^ J dwlnfi A (w) + r ^ +r/^ , 


where k is an integer. For / L (w) > f£ (ui), one can repeat 
similar calculations and obtain the same result. 


(Dl) 


which gives Eq. (l53l) except for a constant tT/2. 

The summation over the discretized counting field, 
Eq. (1571) . can be done for /£"(<*;) < (w) as follows: 


M-l 


£> (ftM+/r(w)e iA *) 
1 /'-/£ 




=Mln/+-^ 


—' n \ f t 

= 1 \ J L 


M-l 

£< 

£=o 


0 i\tn 


=Mln/+ - jr i ( j M(— 


n=l 
. M 


dn (/i (w) M + / L (oj) M ) 


(D2) 
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